Excitation spectrum of a two-component Bose-Einstein condensate in a ring potential 
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A mixture of two distinguishable Bose-Einstein condensates confined in a ring potential has nu- 
merous interesting properties under rotational and solitary-wave excitation. The lowest-energy 
states for a fixed angular momentum coincide with a family of solitary-wave solutions. In the limit 
of weak interactions, exact diagonalization of the many-body Hamiltonian is possible and permits 
evaluation of the complete excitation spectrum of the system. 
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Introduction. Cold atomic gases have a number of ad- 
vantages that make them ideal for the realization and 
study of various physical effects. Topologically nontrivial 
toroidal trapping potentials are of particular experimen- 
tal and theoretical interest, and they have been realized 
by various experimental groups in recent years [l|-l5|- 

Such trapping geometries allow the study of nonlinear 
effects including the excitation of solitary waves and of 
states with circulation. In a classic paper, Lieb consid- 
ered periodic, one-dimensional motion of bosonic parti- 
cles at zero temperature in the limit where the circum- 
ference of the ring, 2'kR, and the number of particles, N , 
tend to infinity with the ratio N / R finite and studied the 
excitation spectrum Q. 

In the present study we consider a mixture of two 
distinguishable Bose-Einstein condensates, confined to a 
torus of a finite radius R, assuming (quasi)-one dimen- 
sional motion. In a previous study [7| we focused on the 
stability of persistent currents as the relative population 
of the two components is varied. Here, we demonstrate 
that this system has a rich variety of interesting proper- 
ties. 

Specifically, we consider solitonic and rotational exci- 
tations of the system. Within the mean-field approxi- 
mation, we evaluate "bound states" of solitary-wave so- 
lutions in the two components. A family of these solu- 
tions actually coincides with the states that result from 
minimizing the energy for fixed angular momentum, the 
"yrast states" in the terminology of nuclear physics [Sj] . 
This identification is general and is thus valid for any 
coupling. This result generalizes that of a recent study 
for a single-component system [9|. The challenge in this 
identification lies in the fact that solitary-wave solutions 
propagate with a fixed velocity while the yrast problem is 
solved for fixed angular momentum. As we demonstrate 
below, for a certain range of angular momentum there are 
"giant" sinusoidal traveling-wave solutions which have a 
constant propagation velocity. We also identify various 
discontinuous phase transitions which occur in the or- 
der parameters of the two components as well as in the 
propagation velocity of the waves. Finally, we go beyond 
the mean-field approximation to describe the exact diag- 
onalization of the many-body Hamiltonian in the limit of 



weak interactions and derive the full excitation spectrum 
analytically. 

Model. The problem that we have in mind is a quasi- 
one-dimensional toroidal potential, with Na and Nb 
number of atoms of species A and B. Within the mean- 
field approximation, the coupled Gross-Pitaevskii equa- 
tions for the order parameters of the two components, 
^A and $B, are 
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where k — {A,B), M is the mass of the two species 
(which is assumed the same for the two components), 
and U = AttU a /{MS) is the matrix element for zero- 
energy elastic atom-atom collisions. (We also assume 
equal and positive scattering lengths for interspecies and 
intraspecies collisions). Further, R is the radius of the 
torus, and 5' is its cross section, with R '3> VS. Finally, 
N = Na+Nb is the total number of atoms, and the order 
parameters are normalized as R J \^k\'^d6 = N^/N . 

Yrast states and solitary-wave solutions. For the yrast 
problem one must minimize the energy of the system 
for some fixed total angular momentum. This can be 
done with the introduction of a Lagrange multiplier. The 
mathematical structure of this problem is closely related 
that for solitary-wave solutions propagating with a con- 
stant velocity. In fact, converting the time derivative into 
spatial derivative, the two problems are seen to be iden- 
tical as we have shown recently in the case of a single 
component [9[. In the present case, the order parame- 
ters of the yrast states, ^ a and 4's, are solutions of the 
following coupled equations: 
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where fl is the frequency of rotation of the density dis- 
tribution of the two components, and l is the operator 
describing the total angular momentum of the system 
(i.e., {La + Lb)Ti divided by the total number of parti- 
cles N = Na + Nb)- Results identical to Eqs. ^ are 
obtained from Eqs. ^ if one makes the solitary-wave 



assumption, i.e., traveling- wave solutions of the form 
$fc(6',i) = e~*'"^-*/'**fc(^ - ^t). The situation is some- 
what more complicated than in the one-component case, 
since there are several families of solutions, corresponding 
to "grey-grey", "grey-bright", and "bright-bright" exci- 
tations. For the repulsive interactions considered here, 
the solitary-wave solutions with the lowest energy, which 
are identical to the yrast solutions, are those for which 
the total density, ua + ?t-b, has the smallest possible vari- 
ation. As we will show explicitly below, these are the 
grey-bright solutions. 

To derive the solitary-wave solutions of the above 
Manakov system [l3] of Eqs. (I2|), we make the ansatz 
ub — KUA + A, where k and A are parameters that can be 
determined from the consistency conditions and the con- 
straints of the problem. The details of this calculation 
will be reported elsewhere. The solutions for the density 
are Jacobi elliptic functions, as is the case for a single- 
component Bose-Einstein condensate confined to a ring 
PEl,i.e., 
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where we write 9 instead oi9 — flt here and below. In the 
above equation /a = sn(a:;|m) and fs — cn(x|TO) are the 
Jacobi elliptic functions, and K{fh) is the elliptic inte- 
gral of the first kind. The phase of the order parameters 
can be obtained from continuity equation, i.e., the equa- 
tions J_ dipk/d9 — 2iTqk, where ipk is the phase of ^fc, 
^fc = \/nk&^'^^i and qk is the winding number for each 
component. 

While the Jacobi solutions are general, they imply 
complicated relationships between the various parame- 
ters. We now focus on several aspects of this problem 
which, although very simple, help to provide a more gen- 
eral picture. In what follows we restrict the value of the 
angular momentum to the range < / < 1. In a single- 
component system Bloch [13| has shown that the more 
general solution and the corresponding energy for any 
other value of I can be evaluated from excitation of the 
center of mass motion. In the present problem with two 
species the same theorem applies Q, and therefore one 
can get the solution for any value of I. 

An exact result. Assume without loss of generality that 
XB = Nb/N is smaller that xa = Na/N. The first 
remarkable result applies when I is in the range < / < 
xb and xa < I l£ i- In this range of I the exact solution 
of this problem is remarkably simple, as it is possible to 
make the total density homogeneous. The easiest way to 
derive this exact solution is to note that when the density 
n^(0) +nB{9) is constant, the coupled, nonlinear Gross- 
Pitaevskii equations, Eqs. ([T]), become linear equations. 
(We stress that identical expressions also follow from the 
solitary-wave solutions of Eqs. ([3]) with k of our ansatz 
equal to —1). More specifically, since the total density is 
constant, the solution of Eqs. ([2]) is 
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where ipm = e""^/v27rR. The corresponding coefficients 
are given by 
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Similar results for the B component are obtained by in- 
terchanging the indices A and B interchanged and with 
dm negative. Here, I is in the interval < I < m. Each 
component has a sinusoidal density distribution. 
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with Ca — cqCto and Cb = dgdm- Since XACoCm + 
XBdodm — 0, the total density UA + nB is indeed constant 
and equal to N/{2TrR). Also 
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where Va — co/cm and Vb — do/dm- The energy per 
particle is given by E/{eN) — 7/2 + ml, where 7/2 = 
2NaR/S is the ratio of the interaction energy of the cloud 
with a homogeneous density N /{2ttR) oi N — Na + Nb 
atoms and the kinetic energy e = fi^ /{2MR^). How- 
ever, since m is an integer, the solution with the lowest 
energy, or equivalently the solution with the lowest veloc- 
ity of propagation u = ili?, is the one with m — \. Thus, 
u = c = h/{2MR), which is constant, simply because 
the dispersion relation scales linearly with the angular 
momentum I, i.e., E/{eN) — 7/2 + I. While the above 
solution is clearly the yrast state, the phase-space con- 
straints do not allow it to have an arbitrary value of / but 
rather restrict it to the range < ^ < xb and xa < ? < 1- 
From Eqs. ([7]) it follows that (for m = 1) the winding 
number of the smaller component B changes at the val- 
ues of the angular momentum li = {1 — -^/l — 2xb)/2 and 
^2 = (1 + VI — 2xb)/2. As we see below, it also changes 
at the symmetry point / = 1/2. (Obviously, it has to 
change in the interval xb < I < xa)- 

Full solution for weak interactions. In addition to the 
exact solution that we have found in the range < I < xb 
and a;^ < Z < 1, the problem also has a simple solution 
in the range xb < Z < 2:^1 in the limit of weak cou- 
pling. Again, the solitary-wave and yrast solutions coin- 
cide. From the point of view of Eqs. ([3]), the parameter 
rh tends to zero in the limit of weak interactions, and 
one can expand the solution in powers of rh. The Jacobi 
functions become sinusoidal since only the single-particle 
orbitals 0o and 0i contribute to the order parameters. In 
this case the parameters of our ansatz are k = —xb/xa 
and A = xbHt^R)- The presence of 4>o and 4>i only is rea- 
sonable if one thinks in terms of the yrast states, where 
minimization of the energy in the limit of weak interac- 
tions is accomplished with these two single-particle states 
of lowest kinetic energy. (Inclusion of (/>_i is also ener- 
getically expensive since it does not have the same sign 
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FIG. 1: The energy of the gas within the mean-field ap- 
proximation (sohd curve) and the lowest eigenvalues of the 
Hamiltonian of Eq. (|13|l (crosses). The dashed curve gives 
the derivative of the mean-field energy, i.e., the velocity of 
propagation of the waves. Here Na = 7 {xa = 0.7), Nb = 3 
{xb ~ 0.3), A'' — 10, and 7 = 2. (Although this is a perturba- 
tive result and 7 ^ 1, we have chosen this value in order to 
make the curves more visible.) The three dots are the points 
with I = h, 1/2, and h- 



as the angular momentum). The order parameters are of 
the same form as Eqs. ([U, vifith 
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The corresponding density and the phase of the two com- 
ponents is given by Eqs. (O and ([7]). Remarkably, in the 
specific case I = 1/2, we find that '^ a — ^/xa{4'o + 
0i)/V2, and *ij == V^(-<^o + 'f)i)/\/2. The point 
/ — 1/2 is singular in the sense that for I — >■ (1/2)^, the 
winding numbers (5^,93) are (0,1) and for I — >■ (1/2)+, 
they are (1,0). The winding numbers thus take the fol- 
lowing values: 



0< ; < (1- VI -2a;s)/2 


{qA,qB) = 


= (0,0) 


(1 - VI - 2a;s)/2 < / < 1/2 


{(iA,qB) = 


= (0,1) 


l/2<;<(l + Vl-2a;s)/2 


{qA,qB) = 


= (1,0) 


(1 + VI - 2a;B)/2 <l<l 


{qA,qB) = 


= (1,1) 



The energy is quadratic in I in this case, 

E/{eN) - 7/2 = ; + -f{xA - m - xb) 



(9) 



(10) 



Its slope, or equivalently the velocity u of propagation, 
varies linearly with I, u/c =1-1- 7(1 — 21). Kt I = xa it 
changes discontinuously by an amount |(5u|/c — ^{xa — 
Xb)- 

While the results presented above for the range xb < 
I < XA are only valid for weak interactions, the more 
general yrast/solitary-wave solutions are given in terms 
of the Jacobi elliptic functions, as described above. For 
stronger couplings, the picture presented for the wind- 
ing numbers remains the same. The density distribution 
may be either sinusoidal or exponentially localized, de- 
pending on the value of the parameters, which determine 



the value of rh. Also as noted above, for m — >■ the 
Jacobi functions become sinusoidal, and for m — > 1 they 
become exponentially localized. 

Exact solution of the many-body Hamiltonian for weak 
interactions. 

What has been presented so far is based on the mean- 
field approximation, which assumes that the many-body 
state is a product state. As we saw earlier, in the limit of 
weak interactions one may work with the single-particle 
states 00 and (j)i only. In this case the many-body Hamil- 
tonian can be diagonalized exactly. We start with the 
Hamiltonian in second quantized form, where Cm, cj„, dm, 
and dl^ with m = 0, 1 are annihilation and creation op- 
erators of the species A and B, respectively. 
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where v = U/{2ttRS). If one considers the algebra of the 
bilinears of the annihilation-creation operators that ap- 
pear in the Hamiltonian, one can recognize two copies of 
the SU(2) algebra and two central elements, correspond- 
ing to the particle numbers of the two species. Let us 
introduce the operators fiAfi = cJcq, ua,! — cjci, ub.o = 
d-odo, fiB.i = ^1^1, and ja — cjco, jb = d[do. The two 
central elements are fiA.Q + fiA.i = Na and fiBfi+fiB,! = 
Nb- If one defines ja,o = {fiA.i - "•A,o)/2, Js^o = ("-s,i - 
hB,o)/2, the two copies of the SU(2) algebra are gener- 
ated by the operators {Ja,Ja,Ja,o}, {Jb,Jb,Jb,o}, where 
[JA,BfijA,B] = JA,B and [Ja,b,Ja,b] = 2JA,B,o- The to- 
tal spin operators (Casimir elements) for the two algebras 

are given by j^^^ = j^bq + (JamJa^b + Ja,bJa,b)/2. 

With the above definitions, the operator that mea- 
sures the angular momentum in units of h takes the form 
L = N/2 + JA,o + JB.o, which allows us to write the 
Hamiltonian in the form 
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Since the possible eigenvalues of ja,Oi Jb.o mn from 
—Na^b/'^ to Na,b/'^, we are in the spin ja,b = Na.b/"^ 
representation of SU(2). Hence j — Na,b{Na,b/2 + 

l)/2. Also, the spin of the j -|-j representation satisfies 
{Na - Nb)/2 < JAB < N/2. This means that the eigen- 
values of the Hamiltonian for a state of definite angular 
momentum L assume the form 



E^eL + v\N\^-l 
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The minimum value of the energy is attained for the mini- 
mum value oi Jab that is compatible with the given value 
of the angular momentum L. One must distinguish three 
cases. For < i < iVs, then jab = N/2 ~ L + k. As a 
result the fcth excited state (with fc = 0, 1, 2, . . . L) is 

Ek-Eo^ {e-v)L + vk{N -2L + k+l). (14) 



For Nb < L < Na, then jab 
thus 



[Na ~-NB)/2 + k, and 



Ek-Eo = eL + v[{Na - L){L - Nb) - Nb] 

+vk{NA-NB + k + l). (15) 

For NA<L<N,we find that jab = L-N/2 + k, which 
implies that 

Ek-Eo = {e + v)L + v[-N + k{2L -N + k + 1)],(16) 

where Eq = vN{N—l)/2 is the energy of the nonrotating 
cloud. The above expressions agree with those of the 
mean-field approximation given above to leading order 
in N ior k = 0. Figure 1 shows the lowest eigenvalues for 
Na = 7, Nb = 3 for the choice 7 = 2. Clearly, the slope 
of the linear parts has a correction which is of order 1/N. 
The case of equal populations Na = Nb = N/2 with 
an angular momentum N/2 is easy to describe and rather 
peculiar. Within the mean-field description, all the an- 
gular momentum is transferred from one component to 
the other as I crosses the value 1/2. Beyond mean-field, 
it is convenient to work in the representation of the Fock 
states |m) = [O^/^"™, 1™) |0™, l^/^-m^^ Here the 

first ket refers to component A and the second to B. 
In this notation N/2 — m atoms of species A reside in the 
state (j)o, etc. It is readily shown that the state 
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is an eigenstate, and actually it is the yrast state, with 
an eigenvalue eN/2 + N{N — 2)z)/2, in agreement with 
Eq. (|15p . This is a highly correlated state and all the four 
eigenvalues of the density matrix are equal to 7V/4. 

Conclusions and overview. A mixture of two distin- 
guishable Bose-Einstein condensed gases which are con- 
fined in a ring potential shows a remarkable collection 



of interesting effects. Although this system is relatively 
complicated, we have succeeded in obtaining a number of 
simple analytic results. In addition we have shown that 
the yrast states coincide with a certain class of solitary- 
wave solutions. 

For a range of angular momentum the two species de- 
velop a quite unusual structure: Although they have a 
long-wavelength and sinusoidal density distribution, ex- 
tending over the whole range of the ring (independent 
of the value of the coupling), they have a finite ampli- 
tude and in that respect differ from the sound waves we 
are familiar with. In addition, their velocity of propa- 
gation is constant, and in the limit of a large ring, this 
velocity tends to zero. When phase-space constraints do 
not allow the angular momentum to take the required 
value, the total density can no longer remain constant 
and the velocity of propagation changes discontinuously. 
The rather complicated solutions in this range of the an- 
gular momentum may be expressed in terms of Jacobi 
elliptic functions. In the limit of weak coupling they take 
again a very simple form. In the same limit we have pro- 
vided the full spectrum of the many-body Hamiltonian. 

The two-component system that we have considered is 
in no sense a trivial extension of the single-component 
system. The combination of two degrees of freedom 
associated with the two species, with a topologically- 
nontrivial motion make this problem very interesting. 
The presence of a second component gives the system 
the freedom to make its total density as homogeneous as 
possible and thus to reduce its interaction energy. This 
is what lies behind many of the peculiar properties de- 
scribed above. The problem examined here is in principle 
more general and thus more complicated than its single- 
component counterpart. In spite of this rich structure, it 
nevertheless has a less complicated solution. We believe 
that experimental confirmation of the effects predicted 
would be both a challenging and a rewarding task. 
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